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ABSTRACT 

We fit model spectral energy distributions to each pixel in 67 nearby {<z>= 0.0057) 
galaxies using broadband photometry from the Sloan Digital Sky Survey and GALEX. 
For each galaxy, we compare the stellar mass derived by summing the mass of each 
pixel to that found from fitting the entire galaxy treated as an unresolved point source. 
We find that, while the pixel-by-pixel and unresolved masses of galaxies with low spe¬ 
cific star formation rates (such as ellipticals and lenticulars) are in rough agreement, 
the unresolved mass estimate for star-forming galaxies is systematically lower then the 
measurement from spatially-resolved photometry. The discrepancy is strongly corre¬ 
lated with sSFR, with the highest sSFRs in our sample having masses underestimated 
by 25% (0.12 dex) when treated as point sources. We found a simple relation to statis¬ 
tically correct mass estimates derived from unresolved broad-band SED fitting to the 
resolved mass estimates: rriresoived = ^unresoZved/(~0-0571og (sS'FR) + 0.34) where 
sSFR is in units of yr“^. We study the effect of varying spatial resolution by degrading 
the image resolution of the largest images and find a sharp decrease in the pixel-by- 
pixel mass estimate at a physical scale of approximately 3 kpc, which is comparable 
to spiral arm widths. The effects we observe are consistent with the “outshining” idea 
which posits that the youngest stellar populations mask more massive, older - and 
thus fainter - stellar populations. Although the presence of strong dust lanes can also 
lead to a drastic difference between resolved and unresolved mass estimates (up to 
45% or 0.3 dex) for any individual galaxy, we found that resolving dust does not af¬ 
fect mass estimates on average. The strong correlation between mass discrepancy and 
sSFR is thus most likely due to the outshining systematic bias. 

Key words: galaxies: fundamental properties, galaxies: statistics, galaxies: stellar 
content 


1 INTRODUCTION 

Exactly how a galaxy’s stellar content evolves over time is 
unknown to us. In general, the build-up of stars—thought 
to be a steady stream of star-formation complemented by 
stochastic bursts triggered by interactions with neighbor¬ 
ing galaxies—is halted at certain times by various quench¬ 
ing mechanisms, such as tidal stripping or feedback from 
e.g. active galactic nuclei. While the overall picture is clear, 
the precise star-format ion history (SFH) of any individual 
galaxy is hidden from us. Over the past decade it has be¬ 
come apparent that a galaxy’s evolution is in large part de¬ 
scribed by its stellar mass. Several tight correlations have 
been found linking the stellar mass of a galaxy to a wide 
range of properties. A galaxy’s brightness and size (Tru- 
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jillo et al. 2006; Williams et al. 2010), rotational velocity or 
velocity dispersion (Faber & Jackson 1976; Tully Sz Fisher 
1977), star formation rate (SFR; Reddy et al. 2006) if on the 
“main sequence” (Daddi et al. 2007; Sawicki 2012b), metal- 
licity (Savaglio et al. 2005; Mannucci et al. 2010), and mean 
stellar age (Gallazzi et al. 2005), are all well correlated with 
the mass of a galaxy’s stellar content. 

Because a galaxy’s stellar mass plays such a vital role 
in our understanding of galaxy evolution, it behooves us 
to understand the estimating methods available as well as 
their limitations. Stellar population synthesis (SPS) mod¬ 
eling (Bruzual Sz Chariot 1993; Bruzual Sz Chariot 2003; 
Maraston 2005; Conroy, Gunn & White 2009; Conroy 2013, 
for a review) is a tool used nearly universally to estimate 
stellar masses in large surveys and at high redshift (e.g. 
Sawicki & Yee 1998; Kauffmann et al. 2003; PerezGonzalez 
et al. 2008; Marchesini et al. 2009; Ilbert et al. 2010; Sorba 
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Sz Sawicki 2010; Maraston et al. 2013). At its most basic, 
an SPS model can be created by combining all the light 
from stellar models in a simple stellar population (SSP)— 
defined as a group of stars all born at the same time, with 
the same metallicity, and with a mass distribution described 
by a chosen initial mass function (IMF)—and tracking how 
the combined light changes as the stars evolve with time. 
Composite stellar population (CSP) models can then be cre¬ 
ated by adding different SSPs together as a function of time 
determined by some hducial SFH. The stellar mass is deter¬ 
mined by fitting a SPS model grid to observed broad-band 
photometry and scaling the stellar mass of the best-fitting 
model. 

The uncertainties in the stellar mass estimate derived 
from broad-band SPS fitting can be large, up to a factor 
of 0.3 in dex when all uncertainties in model assumptions 
are taken into account (Conroy, Gunn & White 2009). This 
large variance does not, of course, account for any assump¬ 
tions that we are not aware of. For example, Maraston et al. 
(2006) found that differing treatment of the thermally pul¬ 
sating asymptotic giant branch (TP-AGB) stage of stellar 
evolution systematically decreased the predicted stellar mass 
of galaxies by a factor of two. Note, however, that there is 
still some controversy over how to correctly treat these stars 
(Kriek et al. 2010). Another model assumption that system¬ 
atically affects the derived stellar mass estimates is one’s 
choice of IMF, since the IMF directly affects the mass-to- 
light (M*/L) ratio of the models. A third systematic bias 
was first noted by Sawicki & Yee (1998) and then empha¬ 
sized by Papovich, Dickinson & Ferguson (2001) who stated 
that star-forming galaxies had the potential to hide underly¬ 
ing older stellar populations by outshining them. The fitting 
procedures preferentially favor matching the large amount of 
flux coming from the younger stellar population, and could 
potentially under-estimate the M*/L ratio by missing the 
relatively low amount of flux emanating from the older stars 
(where a majority of the stellar mass could rest). Unlike the 
two previously mentioned biases, the outshining bias is in¬ 
herent to the broad-band SPS htting procedure, rather than 
emerging explicitly from one’s choice of models. It is thus 
perhaps particularly insidious because it cannot be easily 
corrected for from one study to another. Moreover, because 
the bias would affect star-forming galaxies more strongly 
than quiescent galaxies, its differential nature could act to 
either produce spurious relations or mask true ones. 

The outshining bias essentially arises from our lack of 
knowledge of the SFH of any particular galaxy, and the 
SFH’s influence on the M^/L ratio. If it were possible to 
know the precise SFH to use when creating a model SED, 
then we would know the exact ratio of old to young stars 
and, hence, would not miss any of their mass when fitting 
preferentially to the flux emanating from the younger stel¬ 
lar population. However, since the SFH is unknown, its form 
must be assumed and often parametrized in the model grid 
to account for many possible differing histories. Maraston 
et al. (2010) and Pforr, Maraston & Tonini (2012) used 
mock galaxies to test the effects of mismatched SFHs on 
derived parameters, and found that the stellar masses of 
star-forming galaxies are typically underestimated, in part 
due to outshining, by 0.3 dex for z ~ 2 galaxies and up 
to 0.6 dex for lower redshift galaxies. They found that the 
mass measurement could be better reproduced by using an 


inverted-T SFH (i.e. exponentially increasing), which is pre¬ 
sumably closer to the true SFH of their galaxies. However, 
even if the true SFH is present in the model SED grid, un¬ 
less it is the only SEH present, an outshining bias can still 
occur. Because the flux of younger stellar populations domi¬ 
nates so thoroughly at optical wavelengths, any two models 
with roughly similar young stellar components would have 
roughly an equal chance of being selected as the best-fit 
model, even if one model had the true SEH and the other 
a vastly different one. Outshining can thus be regarded as 
the loss of constraints on the SEH (and thus M*/L ratio) 
because the details are hidden behind bright, young stars. 
Gallazzi & Bell (2009) explicitly examined the uncertain¬ 
ties in M^/L in the case where the true SEH is present in 
the model library, with the assumption that the redshift is 
known and neglecting dust corrections. They were able to 
quantify biases from mismatched SEHs and in the prior dis¬ 
tribution of SEHs compared to the true distribution, which 
they found could be as high as 0.1 dex. 

The extent of any bias from outshining in our obser¬ 
vations of low redshift galaxies has been hinted at, but not 
rigorously determined. Drory, Bender Sz Hopp (2004) com¬ 
pared masses derived from spectra to photometric masses 
and found that for objects with strong Ha equivalent widths 
(i.e. star-forming), the photometric mass was systematically 
lower than the spectroscopic mass by up to 0.15 dex. Zibetti, 
Chariot & Rix (2009) fit SPS models to individual pixels in 
nine SINGS galaxies, and compared the resolved and unre¬ 
solved mass estimates. They also found that the unresolved 
mass estimate (i.e. fitting all the light from the galaxy at 
once as if it were a point source) could be underestimated 
compared to the resolved mass by up to 40% (0.15 dex). 
However, they attribute this effect to unresolved dust lanes, 
rather than star-forming regions. Taylor et al. (2010), on the 
other hand, compared stellar masses with dynamical masses 
and placed an upper limit on any differential bias in stel¬ 
lar mass estimates at less than 0.12 dex. It is clear that a 
systematic bias is present, although the large scatter in any 
individual mass measurement make it difficult to put an ex¬ 
act value on the offset due to outshining. In order to have 
the best understanding of stellar mass estimates, it is crucial 
to better constrain this bias. The large statistical variance is 
stellar mass estimates necessitate the use of a large number 
of objects to penetrate through the noise. 

With this in mind, we follow the approach of Zibetti, 
Chariot & Rix (2009) and compare spatially resolved ver¬ 
sus unresolved mass estimates, as this method only requires 
broad-band photometry which is readily available from a 
myriad of different surveys. We set out to constrain the av¬ 
erage systematic offset from fitting broad band photometry 
of star-forming galaxies to SPS models. We stress that in 
this work we are in no way trying to determine the actual 
amount of stellar material in galaxies, but are instead fo¬ 
cused on a comparative study between the pixel-by-pixel 
and unresolved broad-band SED fitting mass estimates. We 
are solely interested in the effects of outshining as an inher¬ 
ent procedural bias. In Section 2 we describe our data and 
pixel-by-pixel (PXP) SPS fitting procedure, and in Section 
3 we discuss the resulting comparisons between resolved and 
unresolved stellar mass estimates. Throughout this work we 
assume a WMAP7 (Komatsu et al. 2011) flat ACDM cos- 
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mology [Ho = 70.4 km s ^ Mpc Q,m = 0.272) and use 
the AB magnitude system. 

2 METHOD 

2.1 Data 

We used publicly available u,g,r,i, z data from the Sloau 
Digital Sky Survey (SDSS) Data Release 10 (Eiseusteiu et ah 
2011) supplemeuted with publicly available NUV data from 
the GALEX Ultraviolet Atlas (GUA; Gil de Paz et ah 2007), 
which also provides measuremeuts of the spectroscopic red- 
shift, major aud miuor diameter, aud morphological classifi- 
catiou for each galaxy. Although FUV data were also avail¬ 
able from GALEX, we fouud that their sigual-to-uoise (S/N) 
ratio was typically too low to get accurate photometry iu a 
siugle pixel. These six baudpasses should provide more thau 
adequate waveleugth coverage to derive stellar masses from 
SPS httiug, which could be doue usiug ouly the five SDSS fil¬ 
ters^ (e.g. Marastou et al. 2013). The NUV baudpass adds 
a secoud measuremeut below the 4000A break aud provides 
additioual iuformatiou about the youug stellar populatiou 
aud extiuctiou. It acts to break some of the degeueracies 
iuhereut iu the models (see Eigure 1). 

Startiug with the approximately 1000 objects iu the 
GUA which coutaius all galaxies larger thau oue arcmiuute 
iu diameter observed with GALEX, we cross-correlated with 
the SDSS frames usiug the bulk image search fuuctiouality of 
the Sloau website. We theu visually iuspected each match- 
iug frame aud removed auy galaxies that lay too close to 
the border of the SDSS frame. This step eusures that auy 
edge effects, whether from the raw photometry or due to 
couvolutiou (see below), will uot affect a pixel’s hual colors. 
Additioually, we removed auy galaxies that spauued the eu- 
tire SDSS frame siuce there would be uo way to get au accu¬ 
rate backgrouud measuremeut for that frame. The remaiuiug 
galaxies were processed followiug the procedure iu Sectiou 

2.2 aud color images were made from the resultiug aligued 
aud psf-matched images. The color images were agaiu visu¬ 
ally iuspected aud auy galaxies with obvious artifacts (e.g. 
excessive bleediug from foregrouud stars) were discarded. Iu 
total these criteria left us with 67 galaxies iucludiug eight el¬ 
liptical, oue leuticular, aud four irregular galaxies, with the 
remaiuder beiug some form of spiral galaxy. 


^ We compared stellar mass estimates derived using only the five 
SDSS filters with those made with the additional NUV bandpass. 
We found that the inclusion of the NUV had a slight effect on 
mass estimates, increasing them by approximately 7% on average, 
although with large scatter. However, this average difference was 
present in both the resolved and unresolved mass estimates, and 
so does not qualitatively affect our conclusions. We postulate that 
the observed mass difference arises from the NUV limiting the 
amount of unrealistic extinction available to the fit. That is to 
say, we found the pixels (and galaxies) were nearly universally 
fit with higher extinction models when the NUV data was not 
included. The greater amount of extinction for the same amount 
of flux led to greater mass estimates when the NUV data was not 
included. 


2.2 Processing 

Eor each galaxy, we dowuloaded the appropriate SDSS 
frame, mask, aud poiut spread fuuctiou (PSE) which was 
read at the ceuter of the galaxy usiug the read_PSE soft¬ 
ware package provided ou the SDSS website.^ The GALEX 
raw couut, high resolutiou relative respouse, aud flag im¬ 
ages were dowuloaded. We use the GALEX PSE available 
ou the Galex website.^ We created a square cut-out of the 
GALEX images ceutered ou the galaxy with sides equal 
to 1.5 times the galaxy’s major diameter. Although back¬ 
grouud subtracted iuteusity images are provided, we fouud 
through visual iuspectiou of the sky backgrouud images that 
the provided backgrouud maps teuded to overestimate the 
backgrouud of exteuded galaxies, likely due to diffuse radi- 
atiou from the galaxy (see also Zou et al. 2011). We thus 
derived our owu backgrouud estimates iu the raw couut 
cut-out by first maskiug the galaxy aud theu performiug 
a sigma-clippiug fuuctiou iu 50 pixel by 50 pixel square re- 
gious aud httiug a plaue to the backgrouud estimates iu each 
regiou. Siuce couuts from NUV photous are few, the esti¬ 
mated backgrouud value iu each step of the sigma-clippiug 
algorithm was set to be equal to the variauce (rather thau 
the meau) of the values iu each square regiou. This was doue 
to accouut for the lack of symmetry iu the Poissou distri- 
butiou wheu uumber couuts are low; usiug the meau would 
overestimate the backgrouud value slightly. A backgrouud 
subtracted iuteusity map with uuits of couuts per secoud 
per pixel was created usiug the formula 

I=iC-B)/R, (1) 

where C represeuts the raw couut map, B the backgrouud 
map, aud R the relative respouse. We theu created a root- 
meau-square (rms) uoise map for the iuteusity image by 

(7i — ^J + <j‘^l R, ( 2 ) 

where ac is the square root of the couut image (i.e. Pois- 
souiau uucertaiuty) aud cr^ is the staudard deviatiou of the 
backgrouud determiued usiug our sigma-clippiug method 
above. 

The SDSS images are already backgrouud subtracted as 
part of the SDSS pipeliue, but this process over-estimates 
the backgrouud of exteuded objects (Blautou et al. 2005; 
Beruardi et al. 2007). We therefore reiustated the back¬ 
grouud subtracted by the SDSS pipeliue aud performed our 
owu backgrouud estimatiou, similar to what is described 
above, but usiug the meau iu each square regiou as the best 
backgrouud estimate as allowed by the much larger couut 
uumbers. Each SDSS frame had to be couverted to data 
uumbers usiug the provided calibratiou layer of the fits file. 
The rms images were theu made for the SDSS images usiug 
the formula 

ai = (3) 

where I is uow iu uuits of uauomaggies, B is agaiu the back¬ 
grouud map, X is the calibratiou layer value, aud ac, the 
uucertaiuty iu the data uumber image, is uow giveu by 

= VWNUk ( 4 ) 

^ https://www.sdss3.org/drl0/algorithms/read_psf.php 
^ http://www.galex.caltech.edu/researcher/techdoc-ch5.html 
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where C/g is the data number divided by the gain to yield 
the Poissonian variance of the counts, and a dark is the Pois- 
sonian noise from the dark current. 

Having created background subtracted images with cor¬ 
responding noise maps in each of the desired bandpasses, we 
next aligned the SDSS images with the GALEX cut-out and 
resampled the SDSS pixel scale to the GALEX pixel scale 
(1.5 arcsec per pixel) using the SWARP software package 
(Bertin et al. 2002). Although losing spatial resolution seems 
contrary to our purpose, the reasons for degrading the SDSS 
pixel scale are threefold: 1) Artificially chopping up a NUV 
pixel into smaller pixels may not result in realistic spatial 
distribution of the NUV photons (even if interpolation is 
used), which would lead to erroneous colours in each indi¬ 
vidual pixel. 2) Gombining the SDSS pixels together leads to 
a higher S/N ratio which is crucial when striving for accurate 
photometry on a pixel-by-pixel basis. 3) Having fewer pixels 
speeds up computation time of all future steps immensely. 

We also performed a re-alignment and resampling of 
the SDSS PSE by placing the PSE at the same location in 
the frame as the center of the galaxy and using the same 
parameters for SWARP. It is important to have an estimate 
of the PSE prior to resampling, otherwise one is likely to 
have an undersampled PSE and encounter problems when 
trying to estimate it from the resampled image. 

The GALEX PSE was the broadest, so to ensure that 
each pixel contains light from the same physical location it 
was necessary to find a convolution kernel that would trans¬ 
late each image from its current PSE to the broadest PSE. 
We found the transformation kernel by taking the ratio of 
the two PSEs in Eourier space and then inverse-transforming 
the result. We masked any saturated pixels and foreground 
stars (found using the 2MASS catalog; Skrutskie et al. 2006) 
and then convolved each SDSS image with the appropriate 
transformation kernel using scipy’s ndimage convolve rou¬ 
tine with pixel values set to have zero value outside the 
image. The uncertainty was propagated by convolving the 
variance map with the square of the transformation kernel 
to yield the variance map of the convolved SDSS image. 

With the resulting set of background subtracted, 
aligned, resampled, PSE-matched images we then made pho¬ 
tometry catalogs of each pixel with a S/N of at least 5 in 
every SDSS filter. To make a fair comparison between pixel- 
by-pixel SED-fitting result and those using spatially unre¬ 
solved fluxes, we calculated the unresolved magnitudes not 
from aperture photometry, but by summing the flux from 
each pixel included in the pixel-by-pixel catalog and then 
converting to AB magnitude. Doing this is key as it ensures 
that the same amount of light is being compared in each 
case. 


2.3 Models 

We fit the broadband photometry from the six SDSS and 
GALEX bandpasses to model spectral energy distributions 
(SEDs) created using the ESPS population synthesis code 
(Gonroy, Gunn V White 2009; Gonroy V Gunn 2010). Our 
model set is comprised of two component burst models sim¬ 
ilar to those used in Noll et al. (2009) where there is an 
underlying older stellar population combined with a more 
recent burst of star formation. In principle, each of the two 
components (young and old) can have its own SEH, age (i.e. 



Figure 1. Color-color diagrams showing how our model colors 
compare with the observed photometry. The black lines represent 
the 10, 33, 66, and 90% contour levels of the histogram of the 
brightest (Tband) quartile of pixels for all galaxies. The total 
number of pixels binned is 187445. The red points are where our 
models lie in color-color space for a galaxy at 2 ; = 0.0057 which is 
the average redshift of our galaxy sample. The models generally 
show good agreement with the observations, tending to lie on the 
locus in each sub-plot. 


Table 1. Parameters used in two-component model SEDs 


Symbol 

Value(s) 

^old 

0.0016 (~0.084 Zq) 

Zyoung 

0.019 (Z©) 

'^old 

[0.1, 0.25, 0.625, 1.25, 2.5, 5, 10] 

Ryoung 

0.1 

told 

10 

tyoung 

[6, 6.5, 7, 7.5, 7.75, 8, 8.25, 8.5, 8.75, 9, 9.25, 9.5] 

fyoung 

[0, 10-^, 10-3, 0.01, 0.1, 0.25, 0.5, 0.75, 0.9, 0.99] 

E{B - V) 

[0, 0.05, 0.1, 0.15, 0.2, 0.25, 0.3, 0.35, 0.4, 0.5] 


The ages are in units of log years. 


the time since the burst was initiated), IME, and metallicity 
history, as well as various parameters controlling how dust 
absorbs and re-emits light. Here, we describe our chosen 
model parameters, which we found to adequately represent 
our photometric observations (listed in Table 1). 

We parametrized the SEH of the older stellar compo¬ 
nent as a delayed (sometimes referred to as extended) expo¬ 
nential decay (i.e. of the form ~ so that SE at 

first increases linearly with time) with e-folding time Toid 
spaced roughly logarithmically between 0.1 and 10. The 
delayed-T models (Lee et al. 2010) are a good compromise 
between the commonly used r models and the inverted-r 
models of Maraston et al. (2010), particularly for low red- 
shift galaxies which have seen a decline in SER density since 
z ~ 2 (Madau et al. 1996; Sawicki, Lin & Yee 1997). Simha 
et al. (2014) found that delayed-r models were much more 
successful at reproducing the SEH of simulated galaxies than 
plain T models. The SEH of the younger burst is a standard 
exponential decay with Tyoung held constant at 0.1. Each of 
the young and old components also has an age (t). Eor the 
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ngc0628 ngc3367 ngc5719 



Figure 2. Maps showing (from top to bottom) false color ugi, 4000Abreak strength, stellar mass, SFR, and extinction in each pixel for 
three illustrative galaxies from our sample of 67. Areas of low S/N as well as those containing foreground Galactic stars are masked out 
in black. Maps for the remaining galaxies can be found in an online supplement. 
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older component, this was held constant with star formation 
starting 10 billion years ago, while the younger burst had on¬ 
set ages ranging from 10^ to 10® years ago. We combined the 
older and younger components so that the mass fraction of 
the younger component has a value fyoung which we allowed 
to vary between zero and one. The two component models 
used here provide a wide variety of SFHs and are similar in 
construction to those of Noll et al. (2009) which were found 
to adequately represent galaxies in the SINGS survey similar 
to those being studied here. 

We used a Chabrier (2003) IMF for all models and 
set the metallicity of the younger component equal to so¬ 
lar metallicity and that of the older component equal to 
0.0016 (i.e. ~8% Z©). The choice of significantly sub-solar 
metallicity for the older stellar component was motivated 
by the simplistic assumption that the overall metallicity of 
the universe should be lower at earlier times, as well as the 
results of Zou et al. (2011), which showed evidence for a 
double peaked metallicity distribution in NGC 628 with the 
metal poor component peaking at [Fe/H] = —1.15 . Addi¬ 
tionally we found that the inclusion of a sub-solar metallicity 
component yielded a model space that better matched the 
observed pixel-by-pixel colors (see Figure 1). Particularly, 
the g — z color was too red if solar metallicity was used in 
both stellar components. 

For extinction due to dust, we assumed a Galzetti et al. 
(2000) dust law with E{B — V) values ranging between 0 
and 0.5. To save computation time and limit the parameter 
space we arbitrarily enforced models comprised mainly of 
older stars to only have E{B — V) values less than 0.3. 

In total our parameter space is defined by four free pa¬ 
rameters! Told-jtyoung-j fyoung-j and E{B — V). In Figure 1 
we show a comparison between our models at z = 0.0057 
(the average redshift of our galaxy set) and the pixel-by¬ 
pixel photometry for the brightest z-band quartile of pixels 
histogrammed from all galaxies. Although there is a larger 
spread in the observations than the models due to photomet¬ 
ric noise and differing redshifts, our models appear to gener¬ 
ally lie on the locus of the observations and adequately span 
the color-space. It is important to note, however, that since 
we are doing a comparative study of resolved versus unre¬ 
solved fitting, it is not imperative that our choice of models 
actually represent reality to the fullest extent as long as the 
same set of models is used throughout the study. 

In order to determine the best-fitting model to a 
given galaxy’s pixel’s broadband photometry we redshifted 
the model spectrum to the spectroscopic redshift of the 
galaxy and corrected for dust in the Milky Way (Schlegel, 
Finkbeiner & Davis 1998; Schlafly & Finkbeiner 2011). The 
spectrum was then convolved with the broadband filter 
transmission curves to give model magnitudes in each of 
the six bandpasses of interest. We used the SEDfit software 
package (Sawicki 2012a) to generate the model magnitudes 
and then to find the best-fitting model employing a mini¬ 
mization routine. In order to determine uncertainties, SED¬ 
fit performs a number of iterations for each pixel with the 
photometry randomly perturbed by selecting a random vari¬ 
able from a Gaussian distribution with standard deviation 
equal to the uncertainty in the photometry. We performed 
300 iterations and used the median mass of these iterations 
as our final mass estimate. This avoids biases in the best- 
fit estimate (see, for example, Taylor et al. (2011)). 



Figure 3. Ratio of unresolved stellar mass estimate and the PXP 
stellar mass estimate as a function of sSFR. The dashed horizontal 
black line represents a one-to-one correspondence. Spiral galaxies 
(blue points) tend to lie below the one-to-one line, while elliptical 
galaxies (red points) tend to lie on it (on average). The blue line 
shows the average ratio of the spiral galaxies indicating 13 
percent of the stellar mass is missed by fitting all the light as if 
the galaxy were an unresolved point source. The red line shows 
the average ratio for elliptical galaxies, which lies nearly on the 
one-to-one ratio. The shaded blue and red regions show the 1- 
(7 confidence regions for these average ratios as determined from 
Monte Carlo iterations. The purple solid line and shaded region 
show the median and 1-cr confidence region of fitting a broken 
linear relationship to all the galaxies in each of the 300 Monte 
Carlo iterations. The purple dashed line shows a linear fit to the 
data. A clear trend with sSFR is present. Typical error bars are 
shown in the lower left corner. 

3 RESULTS 

In Eigure 2 we show false color images and best-fit model 
property maps for a subset of our galaxies (the full set is 
available as an online supplement). The false color images 
are made using u, g, and z-band images as the blue, green, 
and red colors respectively. The property maps show the 
strength of the D4000 break as a proxy for stellar age (since 
two component models don’t have a well defined age; see 
e.g. Noll et al. 2009), the stellar mass, the star formation 
rate (SER), and the extinction spatial distributions in each 
galaxy. As expected, regions in the spiral arms that are very 
bright in zz-band (i.e. flux generated from younger stars) are 
fit to models with weaker D4000 breaks and SERs that are 
higher than the surrounding population. The mass distri¬ 
bution, however, is relatively smooth, emphasizing that the 
mass-to-light ratio in each pixel is a function of the stel¬ 
lar population, and that older stars comprise the bulk of a 
galaxy’s total stellar mass while not contributing as signifi¬ 
cantly to the galaxy’s total brightness. 

3.1 Unresolved fitting can underestimate stellar 
mass 

In order to determine how well unresolved broadband SED 
fitting can recover the true stellar mass content of a galaxy, 
we calculated the ratio of the unresolved stellar mass divided 
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Figure 4. Same as Figure 3, except showing a false color ugi image for each galaxy. All galaxies have been scaled to be roughly the 
same size, and some of their positions have been slightly altered for clarity. 


by the “true” stellar mass (here assumed to be the sum of the 
stellar mass of the best-fitting model in each pixel). We plot 
this ratio as a function of unresolved specific star formation 
rate (sSFR, SFR per unit stellar mass) in Figure 3. We calcu¬ 
lated the average unresolved-resolved mass ratio for different 
morphologies in each of the 300 Monte Carlo iterations and 
show the median and l-cr spread as horizontal lines in blue 
for spiral galaxies and red for elliptical and lenticular galax¬ 
ies. Doing so reveals that the stellar mass estimates of el¬ 
liptical/lenticular galaxies (red squares) tend to have stellar 
mass estimates that are roughly the same whether resolved 
or unresolved. On the other hand, spiral galaxies (blue cir¬ 
cles) are typically under-estimated by 13% (0.06 dex) on 
average when structure is unresolved. Moreover, there is a 
clear decline in the unresolved-resolved mass ratio for the 
spiral galaxies as sSFR increases. To parametrize this dual 
nature, we fit the broken linear relation 


r hi : X 

I h2 + (hi - h<i)xlp : X > p, 


(which is a horizontal line with intercept hi up to some point 
X — p and has some other slope thereafter) to the galaxies 
in each of the 300 Monte Carlo iterations. The median and 
l-cr confidence intervals of these fits are shown as the purple 
solid line and shaded region in Figure 3. In all Monte Carlo 
instances, there is a clear trend of increasing average mass 
discrepancy with increasing sSFR, with unresolved stellar 
mass estimates underestimated by 25% (0.12 dex) at the 
highest specific star formation rates in our set of galaxies. 
Put another way, to correct the unresolved stellar mass esti¬ 
mate to the pxp stellar mass estimate for these high sSFR, 
one would have to add 33% to the unresolved mass, which 
is much more significant than adding the statistical average 
correction of 15% (0.06 dex). 

Since the break point of the piecewise functions all exist 
at the lowest extremes of our measured sSFRs, separating 
the data into two regimes may be an unnecessary complica¬ 
tion. For simplicity, we fit a simple linear relation to our me¬ 
dian Monte Carlo mass ratios (shown as the purple dashed 
line in Figure 3). This line lies within the shaded l-cr con¬ 
fidence region, and provides a close approximation to the 
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more complex median of the piecewise fit distribution. It 
also allows for a simple formula to convert an unresolved 
stellar mass estimate to pixel-by-pixel one based on the un¬ 
resolved sSFR. The conversion relation is 

^ _ 'kTlunresolved ^ 

“ -0.0571og(s5Fi?) + 0.34’ 
where sSFR is in units of yr~^. 


3.2 Origin of the discrepancy 

In Figure 4 we superimpose false-color images of each of 
the galaxies onto their location in Figure 3. It is apparent 
that the most egregious outliers from the one-to-one rela¬ 
tion (NGC 3190, NGC 4419 and NGG 5719 located in the 
lower left corner with unresolved-resolved mass ratios less 
than 0.65, but also edge-on M98 near the top center) all 
have strong dust lanes with heavy amounts of extinction 
running through them. To try to determine whether resolv¬ 
ing dust or SFR is the dominant driver of the unresolved- 
resolved mass discrepancy, we refit the pixel-by-pixel pho¬ 
tometry but constrained the extinction to be the same as 
determined from the best fit to the unresolved photometry. 
The results are shown in Figure 5 where we plot the ratio be¬ 
tween the fixed extinction PXP mass and the free extinction 
PXP mass versus sSFR. In this figure, the closer the points 
lie to the same mass ratio as Figure 3, the more important 
role dust plays in the unresolved-resolved mass discrepancy, 
whereas points near the one-to-one line are little affected 
by resolving dust. In the cases of NGG3190, NGG4419 and 
NGG5719, the effects of resolving the strong dust lanes ac¬ 
counts for the majority of the difference between unresolved 
and resolved stellar mass estimates, as they again lie well 
below the one-to-one relation. We also find that the large 
majority of the spread in elliptical galaxies and those with 
low sSFR is due to allowing extinction to be a free param¬ 
eter from pixel to pixel. Removing this freedom essentially 
eliminates any difference between resolved and unresolved 
stellar mass estimates for the low-sSFR galaxies. 

What is not present in Figure 5 is any trend with sSFR. 
Even though the mass difference due to resolving dust can 
be drastic for any one galaxy, the deviation appears centered 
on the one-to-one line so that, on average, not resolving dust 
does not systematically bias mass estimates of large numbers 
of galaxies. The strong linear trend seen in Figure 3 must 
then be due to distinguishing the strongly and weakly star¬ 
forming regions and, thus, accounting for the outshining bias 
in broad-band SED fitting. 


3.3 The effect of varying spatial resolution 

In an effort to understand how resolution affects the PXP 
parameter estimations, we artihcially degraded the images 
of the ten largest face-on spiral galaxies by progressively 
doubling the pixel scale. This process was done for those 
ten (of our 67) galaxies with large enough areas that after 
eight degradation steps they were completely unresolved and 
contained within one or two pixels. Eigure 6 shows quartile 
ranges for the ratio of the degraded resolution PXP stel¬ 
lar mass estimate over the full resolution PXP stellar mass 
estimate as a function of physical scale. At the highest res¬ 
olutions, the mass ratio slowly decreases to about 0.97, but 


UJ 1.0 


X 0.8 



log(Free Ext sSFR [yr^ ]) 


Figure 5. As in Figure 3, a mass ratio is plotted versus the sSFR. 
Blue circles are morphologically classified as spirals, red squares 
ellipticals, and green triangle irregular galaxies. Here, the vertical 
axis shows the ratio of the PXP mass found when holding the ex¬ 
tinction value fixed to the amount determined by the best-fitting 
unresolved estimate, to the PXP mass found when extinction is 
a free parameter. 


drops sharply when the pixel scale becomes greater than 2- 
3 kpc/pixel. The mass ratio then stabilizes near the 13% 
under-estimate mark, although with much higher variance 
in the quart ile ranges. It should be made clear that any 
one galaxy does not display a mass ratio that monotonically 
decreases with increasing resolution; the variance of mass 
estimates from SED is too large to show this. It is only by 
looking at an ensemble of galaxies that we can see on aver¬ 
age how the mass estimate behaves as the spatial resolution 
gets poorer and poorer. 

The physical scale of of 3 kpc/resolution element (here 
being pixel because the artihcial degradation makes each 
pixel larger than the PSE) corresponding with the compar¬ 
atively steep drop in mass ratio is much larger than the 
size of typical star forming regions in these galaxies (a few 
hundred parsecs; Efremov 1995; Gusev, Egorov & Sakhi- 
bov 2013). The scale may be more commensurate with the 
width of spiral arms (Kennicutt & Hodge 1982; Seigar & 
James 1998). To adequately separate young and old stel¬ 
lar populations and thus minimize the effect of outshining 
bias, it seems necessary to at least be able to resolve any 
spiral structure in the galaxy. This result may just be an 
artifact of our resolution study sample however, as these 
galaxies all have star forming populations concentrated in 
spiral arms. Note also that several irregular galaxies shown 
in Eigure 4 follow the linear trend due to outshining despite 
the lack of any clear spiral structure. Nevertheless, as spiral 
structure has been observed to have formed much earlier in 
the history of the universe than previously thought possible 
(Law et al. 2012), it is still relevant to know that resolv¬ 
ing spiral arms goes a long way toward mitigating any bias 
due to outshining. Perhaps more generally, one can say that 
pixel-by-pixel SED htting can greatly reduce the effects of 
outshining only when any structure that distinguishes older 
and younger stellar populations is able to be resolved. In the 
case of our low redshift sample, this structure is spiral arms 
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Figure 6. Ratio of degraded resolution PXP stellar mass estimate 
and the highest resolution PXP stellar mass at various pixel scales 
for 10 spiral galaxies with the largest area. The boxes and whiskers 
represent quartiles of the data grouped in bins of width ~ 0.33 
dex in resolution. As in Figure 3, the dashed black line shows the 
one-to-one correspondence and the dashed blue line the average 
mass ratio for spiral galaxies. There is a sharp transition in the 
mass ratio near 3 kpc/pixel. The thin solid coloured lines show 
the individual tracks for each of the 10 galaxies. 

with a physical scale of approximately 3 kpc. If this level of 
resolution is impossible, one can use the conversion relation 
given in Section 3.1. 


4 DISCUSSION 

4.1 Testing Model Space Dependence 

As discussed in Gallazzi & Bell (2009), the model template 
grid used constitutes a prior assumption which can bias 
the mass-to-light ratios determined from broad band col¬ 
ors. In order to test the extent to which our results depend 
on our model space we performed two tests. Firstly, Gal¬ 
lazzi & Bell (2009) found that including a large fraction 
of bursty models could underestimate the mass-to-light ra¬ 
tio of older, smooth-SFH galaxies. Although we found the 
unresolved and PXP masses of ellipticals do not differ on 
average, we nevertheless tested the robustness of our results 
using a smooth SFH. While it is impossible to exhaustively 
test every SFH, we repeated our procedure using r SFH 
models and arrived at qualitatively similar results (namely 
the unresolved mass estimates of galaxies with low sSFRs 
were no different than the resolved mass estimates on av¬ 
erage, whereas star-forming galaxies were underestimated 
by approximately 0.07 dex on average, although without as 
strong a dependence on sSFR). From this we conclude that, 
while our results are not entirely independent of our choice 
of SFH (one could always construct unrealistic SFHs that 
would produce vastly different results, say if all star forming 
models had no underlying older stellar population at ah), 
for any reasonable or commonly adopted SFH our results 
should be qualitatively similar, with unresolved stellar mass 
estimates of star-forming galaxies underestimated by some¬ 
where between 0.05 and 0.08 dex on average. 


Secondly, as examined in Gallazzi & Bell (2009), the 
signal-to-noise ratio of the observations combined with the 
model template prior distribution can lead to biased results. 
In brief, they found that when the model template set is a 
mixture of bursty and smooth SFHs with differing M*/L 
ratios, good data will constrain the hts to only those models 
with the proper M*/L ratio. But, as the amount of noise 
in the data increases, the ht will be more free to explore 
the model space, and hence the determined M*/L ratio will 
depend on the ratio of bursty to smooth SFHs in the model 
template set. Since our stated criteria for a pixel to be in¬ 
cluded in our mass estimation was S/N greater than hve in 
ah SDSS bandpasses, this effect could play a prominent role 
in our results. However, in practice the S/N in g^r^i is much 
better than that in either u or z, and our criteria is essen¬ 
tially equivalent to Su/Nu > 5 (although, in a very small 
percentage of pixels, the z band is the limiting one). The 
large S/N in g,r,i give the hts much less freedom than ini¬ 
tially implied by our selection criteria. To test this, we took 
the unresolved photometry (which has extremely good signal 
to noise) and artihcially increased the uncertainty in each 
bandpass, rehtting the median Monte Garlo mass at each 
step. To realistically approximate the uncertainty growth 
in each bandpass for our data, for each galaxy we created 
histograms of the S/N in each bandpass and found the 
percentile where i is a multiple of ten. We converted this to 
an uncertainty in magnitude for the Monte Garlo hts. Fig¬ 
ure 7 shows the ratio of the increasingly noisy mass estimate 
and the unresolved mass estimate as a function of i^-band 
S/N. The typical r-band S/N ai this stage is also shown for 
reference, and is typically 8-10 times that of the iz-band. The 
binned average (red squares) shows little deviation from a 
one-to-one correspondence. There may be a slight downward 
trend at S/Nu < 10, but it is at the 2% level. Since the low 
S/N pixels are also the least massive, any bias of this level in 
their mass determination has a net effect of much less than 
1% on the total PXP mass of the galaxy. Additionally, since 
the noisy mass would underestimate the true mass, this bias 
would only work to increase the unresolved-resolved mass 
discrepancy found above, if only very slightly. 

4.2 Impact on previously derived relations 

The average offset of 13% for star-forming galaxies is equiva¬ 
lent to approximately 0.06 dex. This is not an outstandingly 
large discrepancy, and any current relationship with stellar 
mass would not change dramatically if corrected for this off¬ 
set. However, the differential nature of the offset (affecting 
higher sSFR galaxies more than those with lower sSFR) can 
compound the effect given the right circumstances. For in¬ 
stance, plotted in Figure 8 in green is the star forming main 
sequence (SFR-M* relation) from Whitaker et al. (2012) and 
isometahicities of the M* — Z — SFR surface from Mannucci 
et al. (2010) in blue. Dashed red lines show corrections to 
these relations based on interpolating the purple line in Fig¬ 
ure 3 of this work. The shallow slope of the Whitaker et al. 
(2012) relation in the SFR-M* log-log plane leads to little 
difference when corrected for outshining. At best one can say 
the slope is a little steeper, requiring slightly faster evolu¬ 
tion from the hatter relation at higher redshift (assuming no 
effect of outshing at higher redshift). However, the effects of 
outshining become more apparent when metahicity depen- 
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Figure 7. Ratio of mass estimates from unresolved photometry 
with artificially increased uncertainties and those with lowest un¬ 
certainties as a function of n-band S/N. The artificially increased 
uncertainty in each bandpass is determined by taking the per¬ 
centile of the S/N histogram where z is a multiple of 10. The typ¬ 
ical r-band S/N determined from these percentiles is also shown 
for reference. The grey points show the mass ratio for each galaxy 
at each noise percentile. The red squares show the average of the 
grey points in bins of 1/10^^ of the data. The black dashed line 
shows the one-to-one correspondence. The mass estimates show 
no bias with increasing S/N within our cut-off limit of 5. 

dence is included in the relationship. The M* — Z — SFR 
relationship exists throughout the SFR-M* log-log plane, 
allowing one to see the larger offsets in high sSFR regions 
(upper left of Figure 8) compared to lower sSFR regions. 
The steeper slope of the isometallicities also lead to out¬ 
shining corrections being larger in the full M* — Z — SFR 
(though still not extremely large). 

The sSFR dependence of the unresolved-resolved mass 
difference can manifest itself in other, possibly less obvious 
ways. For instance, Patton et al. (2013) compare star forma¬ 
tion rates in interacting galaxy pairs with a control sample 
selected to have similar masses as the paired galaxies. Since 
the SFR is enhanced in the paired galaxies, however, so too 
is the sSFR. Thus, a larger correction must be applied to 
the paired galaxy’s unresolved mass estimate than must be 
applied to the control sample, possibly making the control 
sample not representative, ft turns out not to be an issue 
in this particular case, because the maximum SFR enhance¬ 
ment is approximately a factor of 2, meaning the difference 
between the enhanced corrected mass and control corrected 
mass is only 5%, well within their stated control bounds of 
0.1 dex in stellar mass. However, awareness that mass se¬ 
lected samples may not be as homogeneous in mass due to 
differing sSFRs is still important to acknowledge. 

Even within a single galaxy the sSFR-dependant mass 
correction may be important to note. Several works (e.g. 
Allen et al. 2006; Simard et al. 2011; Mendel et al. 2014) de¬ 
compose a galaxy into an inner and outer components (i.e. 
bulge and disc) with SF occuring primarily in the outer re¬ 
gions. The sSFR of the outer region is thus enhanced com¬ 
pared to the inner region, and so a differential correction 
(perhaps up to 0.1 dex) must be applied when comparing the 
masses of the inner and outer components. In general, the 



Figure 8. SFR-M* relation from Whitaker et al. (2012) in green 
and isometallicity contours of the M* — Z — SFR relation from 
Mannucci et al. (2010) in blue. The numbers on the contours are 
in units of 12 + log(0/H). The red dashed lines show corrections 
to these relations based on this work. 

outshining mass correction should be considered for broad¬ 
band SED fitting in any cases with high sSERs, or when 
comparing two sets of data with differing sSERs. 

4.3 Outshining at higher redshifts 

How the systematic bias described above might change with 
redshift is an open question. The sSER of star-forming galax¬ 
ies increases with redshift (Eeulner et al. 2005), and so we 
might expect the outshining effect to become more dras¬ 
tic at earlier times in the universe. Eor example, taken at 
face value, the relation given by the purple line in Eigure 
3 implies that galaxies with sSERs between -8 and -7 yr~^ 
(rates much more common during the era of peak SER den¬ 
sity) would have to be corrected by between 0.1 and 0.2 dex 
respectively. However, as outshining is a product of incor¬ 
rect SEHs, galaxies at higher redshift should be less affected 
because there has been less time for differences in the as¬ 
sumed versus true SEH to accumulate. Put another way, 
there would be fewer older stars (and hence less mass) hid¬ 
ing behind the bright younger stars simply because there 
has not been enough time to form a significant older stel¬ 
lar population. See, for example, Pforr, Maraston & Tonini 
(2012) on the effects of assumed SEH as a function of red¬ 
shift. ft is clear that, at some value of redshift, outshining 
would cease to be an issue as there would not be a significant 
older population of stars to obscure. In the limiting case, the 
first burst of stars formed in a galaxy would essentially be 
a SSP, and thus have no outshining. 

Pixel-by-pixel SED fitting has been performed for high- 
redshift galaxies at z ~ 1 and z ^ 2 by Wuyts et al. (2012) 
using data taken with Wide Eield Camera 3 (WEC3) on 
board the Hubble Space Telescope (HST) as part of the 
Cosmic Assembly Near-infrared Deep Extragalactic Legacy 
Survey (CANDELS; Grogin et al. 2011; Koekemoer et al. 
2011). They found that the unresolved mass estimate from 
these galaxies was robust, and that there was no discrepancy 
between the resolved and unresolved mass estimates. Their 
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PSF matched images had a resolution of 0.18 arcseconds 
FWHM, which corresponds to a physical scale of 1.46 and 
1.54 kpc at redshifts 1 and 2 respectively. From our results 
in Section 3.3 (Figure 6), this resolution should be adequate 
to detect any discrepancy if it existed at high redshift. How¬ 
ever, Wuyts et al. (2012) also further bin pixels to increase 
signal-to-noise, making their effective resolving power much 
worse the the PSF FWHM. Indeed, they state that they do 
find a mass discrepancy of 0.2 dex for z ^ 2 star-forming 
galaxies when no pixel binning is performed, although they 
attribute this to low signal-to-noise pixels having spurious 
mass-to-light ratios. How outshining effects mass estimates 
at high redshift thus remains an open question, one that 
may not be able to be answered until the advent of the next 
generation of telescopes such as the James Webb Space Tele¬ 
scope (JWST) or Thirty Meter Telescope (TMT), which will 
provide roughly 2.5 and 10 times greater resolution than 
Hubble respectively. 


5 CONCLUSIONS 

We fit SPS models to broad-band fluxes in each pixel for 
67 nearby galaxies and compared the stellar mass estimates 
found from summing each pixel’s median Monte Carlo mass 
with the median Monte Carlo mass found when the models 
were fit to all the light from the galaxy at once, as if it were 
an unresolved point source. 

The unresolved mass estimate systematically underesti¬ 
mated the resolved mass estimate, displaying a clear trend of 
increasing discrepancy with increasing sSFR. Galaxies with 
low sSFR {log{sSFR[yr~^]) ~ —12) were consistent with 
no mass discrepancy between resolved and unresolved mass 
estimates, whereas unresolved mass estimates of the highest 
sSFR galaxies (log(sSFR[yr~^]) ~ —8) typically underesti¬ 
mate the mass by 25% (0.12 dex). We provide a conversion 
formula (Equation 6) to correct unresolved mass estimates 
at low redshift based on a galaxy’s sSFR. 

We found the mass discrepancy’s correlation with sSFR 
to be due to outshining, i.e. young stellar populations ob¬ 
scuring older stellar populations behind their bright flux. 
However, several egregious outliers from this relation were 
caused by the presence of strong dust lanes. The presence 
of dust can cause differences between the resolved and un¬ 
resolved mass estimates of up to 45% (0.35 dex). Although 
very significant for individual galaxies, resolving the 2D dis¬ 
tribution of dust did not have any affect on average stellar 
mass estimates. 

By artificially degrading the resolution of the largest 
spiral galaxies, we found that a resolution of at least 3 
kpc is needed to significantly reduce the systematic under¬ 
estimation in stellar mass due to outshining. This scale is 
commensurate with the width of spiral arms in these galax¬ 
ies. 

Although it is unclear how the relation between resolved 
and unresolved stellar mass estimates changes at higher red- 
shifts, we recommend that caution be taken when consider¬ 
ing mass estimates of galaxies with high sSFRs that have 
been derived from broad band SED fitting. Be aware that 
these masses are likely systematically underestimated by 13- 
25% (0.06-0.12 dex) on average due to outshining. The dif¬ 
ferential nature of this effect should be taken into account 


when comparing masses of populations with different sS- 
ERs (e.g. disc versus bulge, “mass-matched” galaxies with 
a range of SERs). In the case of any individual galaxy, note 
that the presence of dust lanes can greatly impact the mass 
estimate (by a factor of ~2 or 0.3 dex), but that the effects 
from dust should average out when considering an ensemble. 
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